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• Sh ! Abstract 
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I. INTRODUCTION 



In radiative processes, the ordering problem of atomic and field operators in the inter- 
action Hamiltonian of bound state QED has been raised since the works by Senitzki [1], 
Ackerhalt et al [2] and others [3]. Behind this discussion is the physical interpretation of 
atomic radiative effects such as the radiative line shifts in spontaneous emission. Alternative 
approaches were proposed in order to elucidate important issues concerning such problem. 
Among them are those based on the complementarity between radiation reaction and vac- 
uum fluctuation effects, which provide a conceptual basis for the physical interpretation of 
different radiative processes. 

In the Dalibard, Dupont-Roc and Cohen-Tannoudji (DDC) formulation, the ordering 
between the operators of the electromagnetic field, considered as a reservoir (TZ), and a 
microscopic atomic system (S) play a fundamental role in the identification of the respective 
contributions due to the reservoir fluctuation (fr) and the self- reaction (sr) [4] - [6]. They 
showed that the symmetric ordering gives a true physical meaning to the (fr) and (sr) rates. 

In this article we study, in the context of the DDC construct, a more general operator 
ordering and its physical significance to a given observable variation rate, more specifically, 
an atomic energy shift. 

We use this analysis to establish a formal connection between DDC approach and a 
closely related treatment, proposed by [6], which is based on the master equation formulation, 
where the physical motivation relies on the classical theory of damping harmonic oscillator 
[7]. Finaly, we discuss the irrelevance of the ordering to an especific interaction Hamiltonian. 

II. THE EFFECTIVE HAMILTONIAN FORMULATION 

In the DDC construct the global Hamiltonian for a coupled system S+TZ is, in the dipole 
approximation, given by 

H^Hs + Hr + V, (1) 
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where Hs is the Hamiltonian of the microscopic system S, Hr the Hamiltonian of the 
reservoir TZ and V the interaction between S and TZ, which we assume to be of the form 
V = —gRS {g is the couphng constant and R and S are, respectively, Hermitian observable 
of 7^ and S) [4] - [5]. 

Following [5] we set that the rate of variation for an arbitrary Hermitian observable G 
of S is given by the Heisenberg equation of motion, and the contribution of the coupling V 
to this rate can be written as 




= -j.[R{t)S{t), G{t)] = gXN{t)R{t) + g{l - X)R{t)N{t), (2) 

coupling 



where N{t) — —{i/h)[S{t),G{t)] is an Hermitian observable of the microscopic system and 
A an arbitrary real number. In the above equation we have used the freedom in the ordering 
of R{t) and N{t), since they commute. 

In order to obtain the contributions of reservoir fluctuation (rf) and self-reaction (sr) we 
perform the following replacement 

X{t)^X\t)+X%t), (3) 

{X — R, S, G) where R^ (resp. and G^) is the solution, to order in g, of the Heisenberg 
equation of motion for R (resp. S and G), corresponding to a free evolution between to and 
t, and R^{t) (resp. and G^) the solution to first order and higher in g. Then, substituting 
(3) in (2) and retaining terms up to second order in g, we obtain 

{t) = -|{(1 - \)R'{t)[S'{tlG'{t)]+\[S'{t),G'{t)]R'{t)} - 

fdt'[S\t'),[S\t),G\t)]]x 

x((l - X)R\t')R\t)) + XR\t)R\t')), (4) 

X ((1 - A)S'(f')[5'(f), G'(()] + AlS'(f), G'(t)]S'(t')). (5) 

Since the rates (4) and (5) contain only free operators, their average value in the reservoir 
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state an gives^ 

{[^fit)) ^ --^^jyC^''\t,t',X)[S\t'),[S\t),G\m, (6) 

X {(1 - X)S\t')[S\t),G\t)] + X[S'{t),G\t)]S'{t')}, (7) 
where we have define g' = \plg and 

A) = \TY^\an{\R'm\^) + (1 - A)i?'(t')i?'(t)}], (8) 
x(^Hi,0 = 'j^^{an{R\^\R\mQ{t-i;). (9) 

The functions C^^^ and x^^^ are statistical functions of the reservoir [8] - [10]. C^^^ is a kind 
of correlation function, describing the "dynamics of fluctuations" of R in the stationary state 
(7ij(io); X*^^^ is the linear susceptibility of the reservoir, determining the linear response of 
the averaged observable {R{t)) when the reservoir is acted upon by a perturbation^. 

The above calculation has shown that the freedom in ordering (2) just reflects in the 
correlation functions given by (8). This result will be explored in section III, where we 
make a connection with the master equation. In order to get a better understanding of the 
aforementioned arbitrariness, we will consider the case of an atomic energy shift. 

In order to find the energy shifts corresponding to the (rf) and (sr) rates we rewrite (6) 
and (7) in a convenient form, namely 

; / \ rf \ (^) 

+ (^) mt,X),[Sit),Git)]] + [Sit),[Y{t,X),Gm)n, 



^Note that the term in the first line of (4) do not contribute to the respective rate since it is linear 
in the absorption and emission operators of the field. 
^In (9) 9 is the Heaviside function, 9{x) = 1 if a; > 0, 9{x) = if a; < 0. 
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=l{[iHeff{t)r,G{t)]U + 



- S(t)[Z"(t,X),G(t)] - [Z'(t,X),G(t)]S(t))n (11) 



where 



{Heff{t)r^'^[Y{t,X),S{t)], (12) 
{HeffiW = ^[Z'{t, X)S{t) + S{t)Z"{t, A)] (13) 

are second order corrections to the Hamiltonian part of S caused by its interaction with the 
reservoir and 

Y{t, X) =J2Qabit){a\S\b) / dr C^^\t, A)e-^--^ (14) 

/oo , , 

drx^^\r)e-'^^^^, (15) 

ao 

dr x(^HT)e-^'^-^ (16) 

with qah = <^a& = (-E'a — E}j)/h and t — t — t'. Following the same point of view of 

[5], expression (12) (resp. (13)) describes the part of the evolution due to reservoir fluctu- 
ations (resp. due to self-reaction) and which can be described by an effective Hamiltonian. 
The second line of expression (10) (resp. (11)) describes the non-Hamiltonian part of the 
evolution of G caused by the reservoir fluctuation (resp. self- reaction) . 



A. The Energy Shifts: Hamiltonian Part 

Corrections (12) and (13) to the Hamiltonian Hs affect S through a shifting in its 
energy eigenstates. Hence, considering a state \a) (which is an eigenstate of Hs) we have 
the following energy shifts 

i5Ej^ = {amffito)r\a), (17) 
{SEj'^^{a\{H,ff{to)r\a). (18) 
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Using expression (12), and noting that 

Y(to)^ C^''\r,X)Sf{to-r)dT, (19) 
Jo 

expression (17) for {SEaY^ becomes 

(SE^f = / ^^""Ht, X)x^''''\r)dr, (20) 

^ J —oo 

where we have introduced a new statistical function, the susceptibihty of the system obser- 
vables 

X^^'^Kr) = ^{a\[Sf{to),Sfito-T)]\a)0{T). (21) 

Prom expression (13) for [HeffY^, we can follow the same steps as those from (17) to (20). 
As a result we obtain 

n'^ r-\-oo 

{6Ear = -V / X^''\r)C^'^''\T, X)dT, (22) 

Z J —oo 

where, again, we have introduced a new statistical function, the "correlation" for the system 
observable 

C(''-){r, A) = ^{a\XSf{to)Sf{to - r) + (1 - X)Sf{to - r)S^{to)\a). (23) 

For future convenience we write (20) and (22) in the frequency space. Using the Parseval's 
theorem we have 

(SEaf = - V / dcv C^^\u;, X)x^''^\oo), (24) 

Z J—oo 

/2 r-\-CO 

{6E,r ---TT dw x^^\u;)C^'''^\u;, A), (25) 

^ J—oo 

where we have used the parity properties of C and x [5] ■ 

Formulas (24) and (25) give us the energy shifts which, a priori, depends on A through 
the "correlation functions", expressions (8) and (23). DDC argued that a true physical 
meaning is obtained by choosing A = 1/2 since then both variation rates, expressions (6) 
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and (7), become hermitian quantities. As a consequence the energy shifts (24) and (25) 
will correspond to the (rf) and (sr) effects. It will be shown in section IV that the above 
assertion becomes meaningless in the case where the interaction Hamiltonian is of the form 

y ~ p • A. 



In this section, we use the previous results to establish a formal connection with the 
approach employed in [6], where the physical interpretation of the energy shifts in term (rf) 
and (sr) are obtained without any reference to operator ordering. In fact, the same energy 
shifts given by (24) and (25) can also be obtained using a matrix approach based on the 
evolution equation for the density operator of the global system S + TZ in the interaction 
picture with respect to Hs + Hr. Hence, following [6], the energy shift for a state \a) oi S 
caused by its interaction with TZ through V is given by 



where is a distribution of probability corresponding to the reservoir average in the sta- 
tionary state Or and |//), \v) are eigenstates of with eigenvalue E^, E^. In (26) V denotes 
the principal value. 

It is directely to see that the matrix element {iJ,,a\V\i', h) in (26) can be factorized in two 
parts, one relative to S and another relative to TZ, 



III. THE MASTER EQUATION APPROACH 




H,i' b 



E^ + Ea — Ej, — Eb 



(26) 




Mi^^ L b 



^ixv + ^ab 



1 



(27) 



In this way, since we know the functional structure of C{ijj) and x((^) [6], namely 



(28) 



X 



(«)(^) = /«)(a;) + zx"(^)(u;), 



(29) 
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V + V- ^ 



(30) 



X"(^)(a;) = I Y.pMR\^)? + - - ^)] , (31) 

and analogous expressions for S (where only pa = 1 is nonzero) , we can make a mathematical 
trick and rewrite the fraction l/(a;^,y + uJab) as 

„ 1 1 



duo X 



X I [V ^- + V ^ 1 + ujab) + (1 - A)5(a; - ujab)\ + 

+ (v^— + V^-)[\5{uj + uj^,) + {l-X)5{uj-uj^,)\\. (32) 

The presence of the A parameter in the correlation functions of equations (24) and (25) 
suggests the above construction. This point is crucial because it clarifies the motivation 
behind the use of (32) in the present context. 

Now, substituting (32) into (27) we obtain: A„ = A^^ + A^^ where 

= -T L W) 

= £ 1^ x'^\u;)C^'-\u;, A). (34) 

Chosing A = 1/2 in the above expression the original formulation given in [6] is recovered 
and the physical meaning of (27) becomes simple and clear in terms of (rf) and (sr) effects. 
This are in complete agreement with the results obtained in the last section, expressions 
(24) and (25). However, in this approach [6] the physical interpretation are borrowed from 
the classical theory of damping harmonic oscillator, without any reference to the operator 
ordering in the correlation functions. In fact, for a class of interaction Hamiltonian of the 
form y ~ p • A the choice of A in (33) and (34) (or equivalently (24) and (25)) becomes 
irrelevant as will be shown in the next section. 
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IV. THE DIPOLE INTERACTION 



In order to show the independence of (33) (or (34)) with respect to A we consider here 
the same case treated in [6] where the interaction Hamiltonian reduces to the expression 



(35) 



where i = x,y, z. Comparing the above expression with the definition ofV {V — —g J2i RiSi) 
we see that g = 1, Ri — Ai{r — 0) (where r = reflects the long wavelenght approximation 
we are assuming) and Si — (epi/m). 

In this context, the final expression for the x-componet of the statistical functions are 

1 



1 



du'uo' 



. 00 + UJ 



UO — 00 



Xr{^) = duo'uo'[5{uo' + uo)- 5{uo' - ^)], 



and 



(36) 
(37) 
(38) 



CTa{^) = E —\{a\p,\h)\M^5{uJab + ^) + (1 - muJab - 



00 



m 



XTaM 



^ 1 ^1 

V + v- 



Xr (^) = E £^l(«KI&)rvr[5(a;„, + a;) - 5(a;„, - a;)], 



(39) 
(40) 
(41) 



where \a) denotes a given eigenstate of Hs- The above expression are obtained following 
reference [6] and tanking into account identity (32). 

In order to obtain the energy shift A^^ we substitute (36) and (40) in (33). Remembering 
that g' = \/2g = \/2, we obtain 

1 



loo 2n [ w? 



V- 



X 



X 



r 1 r'^M ^ 

\ / dJhw'iLniJ)) + 1/2) \\biJ + u;) + (1 - \)b{J - o:)\\ . (42) 

I 377^0^ ^0 J 
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Performing the integral over a;, it is straightforward to see that the contribuction from A and 
that from —A will cancel each other remaning just the contribution from the independent 
term, namely, 

t2 filj7\J ''X 1 



V - + v 



(43) 



where the analogous contribution coming from the y and z components have already been 
added. As we have seen, the above construct clearly shows the independence of (33) with 
respect to A. In the original formulation given in [6] the above expression may indeed 
be divided in two parts: one proportional to the factor 1/2 and associated to the Lamb- 
Retherford energy shift and another proportional to photon nunber and associated to the 
AC Stark effect. 

It must be also noted that the choice A = 1/2 made in the operatorial approach may be 
seen as a fixing parameter which ultimately gives an hermitian character to the contributions 
(6) and (7) allowing a fictions physical meaninging for (24) and (25) in terms of (rf) and 
(sr) effects. However, as remarked by Milonni [11], a true physical meaning is completely 
arbitrary since the observable variation rates are, in fact, unaccessible for the experiment. 



V. CONCLUDING REMARKS 

In this work we have applied to the original formulation of DDC construct a more general 
ordering between the atomic and electromagnetic field (reservoir) operators and calculated 
the energy shift due to the effective Hamiltonian part. The result showed that the freedom in 
ordering expression (2) reflects in the energy shifts (24) and (25) through the A's appearance 
in the correlation functions. 

We have also estabhshed a formal connection between the Effective Hamiltonian ap- 
proach and that based on the master equation theory. This connection was guided by the 
results concerning the general ordering (2) and consequently made explicit the same kind 
of arbitrariness in the physical interpretation of the related variation rates (i.e., the en- 
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ergy shifts (33) and (34)) of section III. In the sequence, it was also discussed that, for an 
interaction Hamiltonian of the form 



the ordering (2) (or (4) and (5)) becomes meaningless and its role on the final result (ex- 
pression (43)) turn out to be irrelevant. As a consequence the Lamb-Retherford shift and 
the AC Stark effect may be obtained without any particular ordering [12]. 

Another interesting result relies on the fact that our procedure still permit us to fix a 
posteriori a suitable ordering which keeps its (rf) and (sr) interpretation, as can be seen by 
looking directly to expressions (33) and (34). However, since the Hermicity of (33) and (34) 
holds only for A = 1/2, such interpretation becomes artificial. 

Once we get a better understanding on the arbitrariness in the operator ordering in DDC 
construct, we expect to find a direct connection with the works by Senitzki, Ackerhalt and 
others. The main idea is to construct a similar structure in Fock space and analyze its 
connection with all possible physical interpretations. 

Another interesting application of the present formalism is a possible generalization of the 
operator ordering in the spirit of q-deformed operator algebras [13], subject of a forthcoming 
work. 



JLT thanks CNPq for partial financial support and the IFT/UNESP for the hospitality. 
LCC is grateful to FAPESP for the financial support. 




(44) 



VI. ACKNOWLEDGEMENTS 



11 



REFERENCES 
[1] I. R. Senitzky, Phys. Rev. Lett. 31 (1973) 955. 

[2] J. R. Ackerhalt, R L. Knight and J. H. Eberly, Phys. Rev. Lett. 30 (1973) 456. 
[3] P. W. Milonni and W. A. Smith, Phys. Rev. A 11 (1975) 814. 

[4] J. Dahbard, J. Dupont-Roc and C. Cohen-Tannoudji, J. de Physique 43 (1982) 1617. 

[5] J. Dahbard, J. Dupont-Roc and C. Cohen-Tannoudji, J. de Physique 45 (1984) 637. 

[6] C. Cohen-Tannoudji, J. Dupont-Roc and G. Grynberg, Atom-Photon Interactions - 
Basic Processes and Apphcations, J. Wiley, NY (1998). 

[7] W. H. Louisell, Quantum Statistical Properties of Radiation, J. Wiley, NY (1973). 

[8] P. Martin, Many Body Physics, Les Houches 1967, edited by C. de Witt and R. Balian, 
Gordon and Breach, NY (1968), p39. 

[9] B. Duplantier, These 3*^ cycle, (1978) Paris (unpublished). 

[10] Expressions of this type can also be found in the context of QED in K. Huang, Phys. 
Rev. 101 (1956) 1173. 

[11] P. W. Milonni, The Quantum Vacuum - An Introduction to Quantum Electrodynamics, 
Academic Press (1994). 

[12] L. C. Costa, master thesis, IFT-D.007/00, IFT-UNESP, Sao Paulo, (2000). 

[13] D. Bonatsos and C. Daskaloyannis, arXiv:nucl-th/9909003 vl (1999). 



12 



